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Abstract 

Within the framework of the Duffin-Kemmer-Petiau (DKP) formalism a more consistent ap¬ 
proach to the derivation of the third order wave equation obtained earlier by M. Nowakowski 
[Phys.Lett.A 244 (1998) 329] on the basis of heuristic considerations is suggested. For this 
purpose an additional algebraic object, the so-called g-commutator (q is a primitive cubic root 
of unity) and a new set of matrices r/ /x instead of the original matrices fi tl of the DKP algebra 
are introduced. It is shown that in terms of these matrices we have succeeded in reducing 
a procedure of the construction of cubic root of the third order wave operator to a few sim¬ 
ple algebraic transformations and to a certain operation of the passage to the limit z —» q, 
where z is some complex deformation parameter entering into the definition of the rj- matrices. 
A corresponding generalization of the result obtained to the case of the interaction with an 
external electromagnetic field introduced through the minimal coupling scheme is carried out 
and a comparison with M. Nowakowski’s result is performed. A detailed analysis of the general 
structure for a solution of the first order differential equation for the wave function z ) is 
performed and it is shown that the solution is singular in the z —>■ q limit. The application to 
the problem of construction within the DKP approach of the path integral representation in 
parasuperspace for the propagator of a massive vector particle in a background gauge field is 
discussed. 
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1 Introduction 


In the paper by Nowakowski [T] devoted to the problem of electromagnetic coupling in the 
Dufkn-Kemmer-Petiau (DKP) theory several rather unusual circumstances relating to a second 
order DKP equation have been pointed out. The first of them is connected with the fact that 
the second order Kennner equation [2] in the presence of an external electromagnetic field is 
only one member of a class of second order equations which, in principle, can be derived from 
the first order DKP equation. Their physical meaning is therefore not entirely clear. Another 
circumstance is connected with the fact that the second order Kennner equation lacks a back- 
transformation which would allow one to obtain solutions of the first order DKP equation from 
solutions of the second order equation. The reason for the latter is that the Klein-Gordon-Fock 
divisor m in the spin-1 case (throughout this work we put h — c—1) 

d{d) = — (□ + m 2 )I + iPud 1 * - — (1.1) 

m m 

ceases to be commuted with the original DKP operator 

L{d) = — ml , (1.2) 

when we introduce the interaction with an external electromagnetic held within the framework 
of the minimal coupling scheme d M —y + ieA i.e. 

[d(D),L(D)}^ 0. 

Here / is the unity matrix; □ = d^d^, d M = d/dx 11 , and the matrices (3^ obey the famous 
trilinear relation 

PuPvP X + P\PuP,x = 9ixuP\ + g\vPn ( 1 - 3 ) 

with the metric g /w = diag(l, — 1, — 1, — 1). Asa result, the analogy of the second order Kennner 
equation to a similar looking Dirac equation is very limited. Whereas in the Dirac case one can 
transform solutions of the second order equation to solutions of the Dirac equation and vice 
versa, such a one-to-one correspondence is not possible in the Kennner case. 

Nowakowski has suggested a way this problem may be circumvented. To achieve the com¬ 
mutativity of the reciprocal operator d(D) and the DKP operator L{D) in the presence of an 
external gauge held we have to give up the requirement that the product of these two operators 
is an operator of the Klein-Gordon-Fock type, i.e. 

d(D)L(D)^-(D 2 + m 2 )I + g{A^ 

where is a functional of the potential A M , which vanishes in the interaction free case. In 

other words it is necessary to introduce into consideration not the second order, but a higher 
order wave equation which would have the same virtue as the second order Dirac equation, i.e. 
a back-transformation to the solutions of the hrst order equation. In the paper [Tj from heuristic 
considerations such a higher (third) order wave equation possessing a necessary property of the 
reversibility was proposed. However, by virtue of the fact that the higher order equation is 
not reduced to the Klein-Gordon-Fock equation in the interaction free case, this leads to the 
delicate question of physical interpretation of the terms in such a higher order equation (it is 
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known that even for the second order Kennner equation there exists this kind of problem). In 
particular, this is concerned with the interpretation of parameter m as the mass of a particle 
since this is the only possible interpretation in the free case when the following equality, 

d(d)L(d) = -(D + m 2 )/ (1.4) 


holds. 

It should be noted that the divisors with a minimal electromagnetic coupling and, in par¬ 
ticular, for the spin-1 case 

d(D) = — (D 2 + m 2 )/ + i/TD^- —/y^DGD 1 ' (1.5) 

m m 

were first introduced into consideration in the earlier papers by Nagpal [5j, Cox )bj, and Krajcik 
and Nieto [7jj. The divisors have been intensively used in analysis of causality violation in higher 
spin theories in the presence of an electromagnetic field. The suggested divisors represent merely 
a straightforward generalization of the well-known operators of Takahashi and Umezawa HU by 
the replacement d M —> D^(A). However, a question of commutativity of the generalized divisors 
with the initial first order operators L(D) in these papers was not discussed at all, although 
this can be of certain importance. Further, in the papers mentioned above the questions of 
causality were discussed on the basis of analysis of a product of two operators d(D) and L(D). 
In particular, in the spin-1 case when we take the divisor in the form (11.51) . in the product 
d(D)L(D) the principle part of interacting and free wave operators remains the same as it 
was defined by Eq. (11.4j) . This is connected with the fact that the terms of the third order 
in derivatives reduce to the terms of the first order by using the trilinear relation (11.31) and 
thus the effect of electromagnetic interactions (or nonderivative coupling) occurs only in lower 
derivatives. The resulting field equation remains equivalent to a hyperbolic system with light 
cone as ray cone, the same holds in the interacting and free cases. Therefore, it is concluded 
that the spin-1 field even in the presence of electromagnetic field in the system possesses only 
causal modes of propagation. 

The situation can qualitatively change if as d(D) one takes a divisor such that 

[d(D),L(D)} = 0, 

for example, the divisor suggested by Nowakovski [I]. In this case in the product d(D)L(D) 
in accordance with formula (16.101) in Section 7, the principle part of the interacting and free 
wave operators will be already the third order in derivatives, instead of (11.41) . and the terms 
with the nonderivative coupling remain the terms of the first order. The question of whether a 
change of the order of the principle part of wave operator leads to a change of the propagation 
properties of the equation 

d(D)L(D)ip(x) = 0, 

generally speaking, has to be the subject of separate research. 

Further, in addition to the absence of required one-to-one correspondence between solu¬ 
tions of the second order Kennner equation and the DKP equation, one can point out one 
more negative consequence of the noncommutativity of the divisor (11.51) with the operator 
L(D ) = — ml. The lack of commutativity does not give a possibility within the frame¬ 

work of the DKP theory to construct the path integral representation for the Green’s function 
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of a spin-1 particle in a background gauge field in a spirit of the approaches developed for a 
spin-1/2 particle (see for example Els]). Having obtained all the necessary expressions, this 
very interesting question will be discussed in more detail in section 8. 

The purpose of this paper is to give a systematic way of deriving the third order wave 
equation within the framework of the massive Duffin-Kemmer-Petiau theory in the free and 
interacting cases. However, first of all it should be noted that the wave equations of the third 
order in derivatives, as applied to the problems of classical and quantum held theories, for any 
length of time have drawn attention of researchers for various reasons. Below, we give a number 
of examples related somehow to our problem. 

In the papers by Finkelstein et al. ra in constructing the theory of the nine-dimensional 
ternary hyperspin manifold the so-called trine-Gordon equation, the unique scalar wave equa¬ 
tion of least differential orderQ 

[det(<9) — im 3 ](p = 0 

was suggested. Here, det(<9) is the determinant of a 3 x 3 matrix composed of partial derivatives 
in coordinates. The authors have also performed an analysis of the corresponding dispersion 
relation for plane waves and have suggested the generalization to the case of a minimal interac¬ 
tion with a gauge held. The questions close to this research were considered independently in 
the papers by Solov’yov et al. HU, where the general algebraic theory of the Finslerian spinors 
was constructed. The generalized Dufhn-Kemmer-Petiau equation for a Finslerian 3-spinor 
wave function of a free particle in the momentum representation was also suggested there and 
it was shown that each of these 3-spinor components of the wave function (i r , /3$) satishes the 
Finslerian analog of the Klcin-Gordon-Fock equation 

(G A bcP A P B P C ~ M 3 )i r = 0, r = 1, 2, 3 


and a similar equation holds for the /h, components. Here, Gabc is a symmetric covariant 
tensor of the third order rank that plays a role of the metric tensor in the nine-dimensional 
linear Finslerian space. 

Further, in the works by Yamaleev HM3 an attempt has been made to construct in a 
systematic way the foundations of quantum mechanics on cubic forms (or even more generally, 
polylinear forms). The mathematical basis of the construction would become the cyclic algebras 
of N > 2 degree. The cyclic algebra with respect to cubic forms here plays a role like the 
Clifford algebra with respect to the quadratic forms. In particular, the cubic generalization of 
the standard relativistic relation between energy £, momentum p = (p \, P 2 PP 3 ) and mass m was 
suggested in the following form: 

3 

(£ - qm){£ - q 2 m)(S - m) = '^pf - 3 pip 2 p 3 , (1.6) 

2 = 1 


where q is a primitive cubic root of unity 


q = e 2ni/3 = 
q 2 = e 4 ™ /3 



(1.7) 


X A11 formulas cited below, up to Eq. m , are given in the notations of the authors of the corresponding 


works. 
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As an analog of the Klein-Gordon-Fock equation (for any of three possible correspondences: 
£ —» 9kd/dt,pi —$■ Qkd/dxi , where Qy. = (g, g 2 ,1), k = 1, 2, 3 in Eq. (11.61) 1 the following equation 
of the third order 



d s 


<p — m 3 p = 0 


dxidx 2 dx 3 

was used. The author also suggested cubic analogs of the Dirac equation and the nonrelativistic 
Schrodinger equation. 

The third order wave equations arise, however, not only in the generalization of quantum 
mechanics and quantum held theory to more abstract spaces in the foundation of which not the 
quadratic forms of various type (an interval, relativistic relationship between mass, energy and 
momentum and so on), but the forms of one degree higher are laid. These equations arise also 
within the framework of generally accepted physical theories for solving quite concrete problems. 
Thus, one of the first generalizations of this kind in the context of quantum electrodynamics can 
be found in the paper by Pais and Uhlenbeck US- The latter have considered the generalization 
of the Dirac equation to the multimass Dirac equation like 


N 

JJ ('y^d tt + m j )'i/j{x) = 0 , 

3 = 1 


where N is an integer, which in particular can take the value 3. The purpose of this general¬ 
ization of the spinor held equation to equation of higher order was to eliminate by this means 
the divergent features in quantum held theory. 

Further, in the paper by Barut et al. [IB] another version of the generalized Dirac equation 
of the third order in derivatives describing particles with spin 1/2 and three mass states was 
suggested. The aim of this paper was to treat in a unified fashion all known at that time 
leptonic matter. The case when one of the states is massless (neutrino) and the corresponding 
generalized equation takes the form 


\ioti {7 • d ) + a 2 D — * 0 : 3(7 • d)\3\ijj{x) = 0 


was studied in more detail. The parameters a±, a 2 , and 03 are related to the electron mass m e 
and muon one through the relations 01/0:3 = m e m^, a 2/0:3 = m e + It is interesting to 
observe that the term of the third order in derivatives in the equation above has the structure 
similar to the corresponding term on the right-hand side of Eq. (II.lip for the spin-1 case, which 
will be discussed further. In the paper by Kruglov m this model was analyzed for the case 
when all three states are massive. As was noted in the last paper, such higher order differential 
equations may be treated as effective equations and represent a nonperturbative approach to 
quantum held theory. 

In the spin-1 theories by Joos [18], Weinberg [IB] and Shay et al. [20] it was shown that 
the corresponding wave function i^{x) in the interaction free case must satisfy component by 
component not only the second order Klein-Gordon-Fock equation, but also the third order 
wave equation of the type 


^liv A 


d 3 ip(x) 

dxndx u dx\ 


m 3 ijj(x) = 0 . 


Here the 8 x 8 matrices 7 ^ are defined in terms of the spin-| matrices s t , i = 1, 2, 3 and obey 
algebraic relation representing the spinA generalization of relation for the Dirac matrices. 
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Finally, we can also mention that in the familiar formulation of Bliabha I 2 U (see also [ 22 ]) 
of the multimass high-spin theory, for the spin-1 case we have instead of the Klein-Gordon-Fock 
equation the third order one 

(a ■(?)(□ — m 2 )ip = 0 . 

The extra differential factor (a • d ) conies from the subsidiary components. 

The examples given above show that the higher order systems, in particular the third order 
ones, might themselves have some applicability in field theories. 

Before proceeding with the formal development of the construction of the third order wave 
equation within the framework of the DKP approach in the presence of an electromagnetic 
field, it is necessary to ask what form the equation should have in the interaction free case. It 
is necessary to have at hand a certain simple rule of deriving this equation (and perhaps, the 
equations of higher order for high-spin cases). Here, we attempt to follow as close as possible 
the free Dirac theory added by some considerations of algebraic character. For this purpose 
let us introduce a set of the square roots of unity: (A, 1), where A = —1 is the primitive square 
root. Then it is obvious that 

(iy^cP — A ml) [i^ v d v — ml) = — (□ + m 2 )I. (1.8) 

Let us state a question of defining such a matrix O that the following equality holds 

\0{i y^ck 4 — ml)] [(D(i r y„d l/ — ml)] = — (□ + m 2 )I. (1.9) 

In fact it represents a solution of the problem of constructing the square root of the Klein- 
Gordon-Fock operator. Here, the answer is known. As such a matrix one has to take 

O = ±«y 5 , 


where 

y5 = ^ ^ivTivixi*, il = i- (i-io) 

Thus we can consider that the expression on the left-hand side of (11.81) gives us the rule for the 
determination of the right form of the second order wave operator (the right-hand side of (I1.8P ) 
and in turn the expression on the left-hand side of Eq. (II .9j) gives its square root and thereby 
the problem is reduced to the construction of an algorithm of calculating the matrix O. 

If one consider as a guiding principle the considerations above, then the next step will be 
the following extension: as a basis we take the cubic roots of unity ( q , q 2 , 1 ), where the primitive 
roots q and q 2 are given by the formulas (11.71) . and as the spin matrices we take the /3-matrices 
of the DKP algebra. It is an easy matter to verify that an analog of equation (11.81) will be that 
in the following form (cp. with (11.61) '): 

(i/3 M cP — qml) (if3 u d v — q 2 ml) ( i/3\d x — ml) = — i □/3 /i <9 /i — vrC’I. (1-11) 

On the right-hand side of (11.111) we now have the differential operator of the third order, 
which we take as a “genuine” expression for the third order wave operator. It is precisely this 
expression that arises from the Nowakowski’s third order wave equation [T| in the limit when 
we switch off an external electromagnetic field. In deriving (II.lip one of the properties of the 
roots of unity, namely 

1 + g + g 2 = 0 (1.12) 
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and the identity 


^M x d^d x = □ 0^ (1.13) 

valid in view of the algebra of the /3-matrices, Eq. (11.31) , were used. Besides, we have taken into 
account the fact that the mass term ml is diagonal and commutes with everything. We note 
that such an approach was used in the papers by Kerner [23] devoted to a generalization of 
supersymmetry based on z/ 3 -qraded algebras, more exactly in the construction of the operators 
whose trilinear combinations yield the supersymmetric generators (cubic root of the supersym¬ 
metry (SUSY) translations). 

Further we can state a question of defining a matrix A such that the following relation holds: 

[^(i/^cU — ml)] [A(if3 u d 1 ' — ml)] [ A(i/3\d x — ml)] — — i — 0/3^8^ — m 2 I. (1-14) 

The latter solves the problem of calculating the cubic root of the third order wave operator. 
In this paper we have attempted to answer this question by using a very rich apparatus of the 
matrix algebra in the DKP theory added by new structures generated by algebra of the cubic 
roots of unity. We have also performed a generalization of the resulting equations to the case 
of the presence in the system of an external electromagnetic field. 

The paper is organized as follows. In Section 2 the construction of cubic root of the second 
order Klein-Gordon-Fock operator within the framework of DKP formalism is considered. This 
problem has a purely auxiliary character. However, a number of expressions derived here are 
of decisive importance for subsequent research. Section 3 is devoted to the construction of the 
cubic root of the third order wave operator. For this purpose an additional algebraic object, the 
g-commutator representing a generalization of the usual commutator by entering a primitive 
cubic root of unity q into the initial definition, is introduced. This new algebraic object has 
allowed us to remove not only the terms linear in derivatives, but also the quadratic terms as it 
is required by virtue of the definition of the third order wave operator, Eq. (11.141) . However, at 
the same time it is found that the necessary term cubic in derivatives, on the symmetrization, 
vanishes identically. 

In Section 4 a new set of matrices r]^ instead of the original matrices /3 M is introduced. It is 
shown that these matrices possess rather nontrivial commutation relations which enable us to 
reduce the problem of the construction of the desired cubic root to a number of simple algebraic 
operations. On the basis of these matrices the reason for vanishing the term of the third order in 
derivatives is analyzed and a way to overcome this problem is suggested. Section 5 is concerned 
with the discussion of various properties of the 77 -matrices: commutation relations, the trilinear 
relation (the analog of the trilinear relation for the /3-matrices), the behavior on Hermitian 
conjugation, etc. In Section 6 an extension of the results of the previous sections to the case 
of the presence in the system of an external electromagnetic held is performed. The detailed 
comparison of the expression for the third order wave operator with a similar expression earlier 
obtained by Nowakowski [I] is given. In Section 7 an analysis of the general structure for a 
solution of the Erst order differential equation for the wave function ^(x^z), where z is the 
deformation parameter is performed. It is shown that this solution is the singular one in the 
limit z —y q. In Section 8 a question of a possible application of the results obtained to the 
problem of the construction within the framework of the DufEn-Kennner-Petiau formalism of 
the path integral representation for the propagator of a vector particle in a background gauge 
held is discussed. In the concluding Section 9 the key points of our work are specified and the 
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massless limit of the third order wave operator is briefly discussed. 

In Appendix [A] all of the necessary formulas of the DKP algebra of the matrices b ;/ are 
listed. In Appendix [B] a procedure of the construction of a certain matrix A is presented. This 
matrix formally can be considered as a matrix analog of the primitive cubic roots of unity, i.e. 
a set of the matrices {A, A 2 , A'*’ = ^/) satisfies the properties identical to those for a set of 
the cubic roots of unity: (g,? 2 ,? 3 = 1). In Appendix [Cl a complete proof of vanishing cube 
of matrix differential operator, where matrix is defined through the deformed commutator is 
produced. In Appendix [D] the details of the proof of trilinear relation of the type (11.31) for a 
new set of the matrices are given. Finally, in Appendix [E] the proof of the identity (16.31) for 
a product of three covariant derivatives is presented. 


2 Cubic root of the Klein-Gordon-Fock equation 

Before proceeding to the problem stated in Introduction, we first to consider a question of the 
construction of cubic root of the second order massive Klein-Gordon-Fock operator. This prob¬ 
lem in the general statement has been investigated by Plyushchay and Rausch de Traubenberg 
in the paper m- Here, we examine it again and look how far we can proceed in solving this 
problem while remaining within the framework of DKP formalism only. 

Let us now turn to Eq. (II. 14)1 . but instead of the third order operator we put the Klein- 
Gordon-Fock operator on the right-hand side 

— ml)] [A(i/3 l/ d L ' — ml)] [A(i/3\d x — ml)] = —(D + m 2 )/. (2.1) 

One can somewhat simplify the problem if one takes the operator on the right-hand side in the 
factorized form 

— (□ + m 2 )I = d(d) ( M — ml ). 

We recall that 

d(d) = ml + i/3 ^ + (2 - {/^, A,}) - 

is the Klein-Gordon-Fock divisor in the spin-1 case; { , } designates anticommutator. By virtue 
of this factorization we can examine instead of (12.11) the following equation: 

A[if3^d^ — ml)A[i/3 u d l/ — ml) A = d(d ). (2.2) 

By equating the coefficients of partial derivatives we obtain a system of algebraic equations for 
the unknown matrix A: 

A 3 = — I, (2.3) 

m 

Af3 ll A 2 + A 2 (3 ll A = --f3 ll , (2.4) 

m 

A^A(3„A + APvAPpA = - ^ [2 g^I - {/9„, &}]. (2.5) 

These equations can be paired with the corresponding equations in the paper [24], if we correlate 
the generators g M and g introduced in [24] with the matrices A and (3^ by the rules 


g,i ~ mA(3n, g ~ mA. 







In this case, Eq. (12.31) (up to a sign) will correspond to the first equation of the system (12.31) in 
the paper Eza ,Eq. (12.41) will correspond to the second equation of the same system (or Eq. (12.6ft ). 
and (12.51) corresponds to the third equation. 

Before turning to solving the matrix equations (12.31) — (12.5j) , we make a few comments of a 
general character. Equations (12.3ft and (12.4ft are universal in the determinate sense. The former 
defines the mass term on the right-hand side of the equality (12.11) (or (11.141) 1. and the latter 
enables us to get rid of the term of the first order in the derivatives in (12.11) (or (11.14ft ). The 
universality of these matrix equations lies in the fact that they must be satisfied in any case 
irrespective of that we take as the right part: either the right-hand side of ( 12 . 1 1 ) or the right- 
hand side of (11.141) . As will be shown below, Eqs. (12.3ft and (12.41) uniquely define the required 
matrix A (more exactly, to within the choice of one of three roots of the cubic equation for some 
parameter a, see Eq. (I2.8[) below). An explicit form of the matrix A and also the equalities 
(12.31) and (12.4ft to which it satisfies, are of fundamental importance for further presentation. 
The third equation (12.51) is not already universal and completely depends on the specific choice 
of the right-hand side in the equalities of the (12.11) type. This equation must be identically 
satisfied. If not, we come to the contradiction. 

Let us now introduce the matrix oj setting by definition 

w = \ (2-6) 

This matrix plays an important part in further consideration. It was introduced into DKP 
theory for the first time by E. Schrodinger [25]. Here, we follow the notation used in the 
works by Harish-Chandra [26], where the properties of the c o matrix were studied in detail. In 
Appendix A we give all necessary relations for the algebra. Let us note only that the 
matrix c o is identically zero for the spin 0 (five-dimensional irreducible representation of the 
DKP algebra). Therefore, only the ten-row representation needs to be considered. 

In spite of a formal similarity between definitions 75 and oj matrices (it is worthy of special 
emphasis that in the latter case the factor 1/4 stands rather than 1/4!), Eqs. (11.10ft and (12.61) . 
as the matrix A we cannot simply take Really, for example, on the left-hand side of 

Eq. (12.31) in view of (lA.lft we will have A a ; 3 = A u 7 ^ A / 

We seek the matrix A in the form of the most general expansion in powers of uj: 

A = al + f3u) + 7 cn 2 , 

where a, /3, and 7 are unknown, generally speaking complex, scalar constants. By virtue of the 
property (lA.ll) it is easy to find that 

A 2 = a 2 1 + (2 af3 + 2/3 7 )u; + (/3 2 + 7 2 + 2 a^)oj 2 (2.7) 


and further 

A 3 = a 3 1 

+ [a(2 otf3 + 2 /^ 7 ) A ot 2 /3 + 7 (2a/3 + 2 /^ 7 ) + (/3 2 + r y 2 + 20 : 7 )] w 
+ [ a(/3 2 + 7 2 + 20 : 7 ) + /3(2a/3 + 2 /L 7 ) + a 2 7 + 7 (/ 3 2 + 7 2 + 20 : 7 ) ]co 2 

= —I. 

m 
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The foregoing expression enables us to reduce Eq. (12.31) to a system of three algebraic equations 
for unknown scalar constants, the first of which defines the parameter a: 


a 


3 


1 

m 


( 2 . 8 ) 


Two other equations follow from vanishing the expressions in square brackets. However, instead 
of these equations it is convenient to consider their sum and difference, which after simple 
algebraic transformations can be recast in a more convenient form 


(/3 + q) [(/3 + q) 2 + 3a(/3 + q) + 3a 2 ] — 0, 

(P ~ l) [(P ~ if ~ 3 a(P - q) + 3a 2 ] = 0. 

As solutions of these equations we take solutions of the quadratic equations in square brackets 
for the variables (/3 + q) and (P — q), namely 


{P + i)± = (“ ± 

(/3-7) ± =(5±»^)a. 


We return to the obtained solutions just below, and now we pass on the second matrix equation 
(12.41) . By using the properties (IA.1D - (IA.2D . we get 


A^A 2 + A 2 f3^A = 

— [ 2 a 3 + a 2 q + a(/3 2 + q 2 + 2 aq)]/?^ + [a 2 /? + a(2a/3 + 2/3q)j (tn/3^ + u) = — -— /3^, 
which due to (12.81) gives us the second system of algebraic equations: 

a(/3 2 + q 2 + 3 aq) =-, 

m 

a(3 a(3 + 2/3q) = 0. 


Since a^O and considering f3 ^ 0, from the very last equation we obtain 



Taking into account this fact, from the first equation we derive /3 2 = —(3/4) a 2 or P = 
±i(\/3/2)a. The solution obtained for the parameter q is not in contradiction with the solu¬ 
tions m- As the f3 parameter one can take either of the two solutions ±i(\/3/2)a. It is also 
consistent with the solutions (12.91) . For definiteness let us fix the sign +, i.e., we set 


o . V3 
p = i — a. 
2 


The matrix A in this case takes the following form: 


A = a 




( 2 . 10 ) 
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An explicit form of the matrix A 2 can be obtained by squaring (12.1011 or by making use of (12.7|) . 
Here we have 

1 ( 2 - n ) 

Note that A 2 is different from A not only by the extra dimension factor a, but also by the 
opposite sign before the second term (in fact, here we have the second possible value of the pa¬ 
rameter j3). An explicit form of the matrices A and A 2 hints that they are mutually conjugated. 
Indeed, if one takes into account the Hermitian character of the cu-matrix: 



U!' = uj , 


then the following two relations are true 

= m 1/3 A 2 , AA f = A f A = I. (2.12) 

In principle, one can avoid the mass multipliers if one overdetermines the matrix A, but we do 
not do it. 

Furthermore, we can draw an interesting parallel between a set of matrices ( A , A 2 , A 3 = A /) 
and a set of cubic roots of unity: (g, g 2 ,g 3 = 1), Eq. (11.71b In the latter case the following 
relations, which are similar to (j2. 12f) . hold, 

q =q , qq = 1, 


where q* is the complex conjugate of q. However, the cubic roots (g, g 2 ,1) possess one more 
important property (11.121) . whereas for the matrix set (A, A 2 , A /) we have 

I + -A + ^A 2 = 3(I-cu 2 ). 
a cr 

Nevertheless, it is possible to redefine the matrix A —\ A such that the following equality will 
be held 

I + AA + -A-A 2 = 0 

a cr 

and at the same time the properties 

A 3 = — I, A ] = m l ' 3 A 2 

m 

will be survived. Since the matrix A will not play any role later, we give its explicit form in 
Appendix B. 

Now we turn to analysis of the remaining equation (12. 5|) . By making use of an explicit form 
of the matrix A, Eq. (I2.10K . and the properties (1A.4D - (1A.7D . we obtain for the left-hand side of 

m 


AfinAfivA + AfivAfi^A — — 


m 


1 \/3 3 3 

2 A i - 9fj .A — 9 in ,uj — {/^^, /3[y}cu 


(2.13) 


Comparing the expression with the right-hand side of (12.51) . we see that their matrix structure 
is sufficiently close to each other. The main difference with (12.51) is the presence of the term in 


11 
























(1 2.13 p linear in the matrix uj. We can remove this term if we slightly complicate the left-hand 
side of the initial expression ( 12 . 2 [) . namely, we present it in the following form: 


1 (AL(8)AL(8)A + mA 2 L(8)A 2 L(8)A 2 ), 


(2.14) 


where the operator L(d ) was specihed by Eq. ( 11 . 21 ) . It is not difficult to see that the first two 
equations ( 12 .3|) and (12.41) remain unchanged and instead of (12 . 51) now we have 

— [{A(3 ll A(3 u A J rA(3 l/ A(3 fl A) J rm(A 2 f3 fl A 2 (3 l ,A 2 - J rA 2 (3 l/ A' J fi fl A 2 )^ = — — (2$^—{/3^,A,}). (2.15) 

By using an explicit form of the matrix A 2 one can see that the expression on the left-hand 
side of (12.151) completely coincides with (I2.13P except for cancellation of the term linear in oo. 
Thus, the matrix equation (12.151) leads to fulfilment of the following equality: 

2 {An “b 2 9nv^ ~ 2 m = 2 ~ {(3[n/3 U }■ 

The relation is inconsistent. To verify this, it is sufficient to contract it with g for example. 
With the relations (1A.9|) we result in a contradiction 



One can look at the problem in a different way. As we know in the interaction free case the 
divisor d(d) commutes with the operator L(d ), i.e. 


[■ d(d),L(d )] = 0 . 

Let us substitute now the operator (12.141) instead of the divisor d(d). The result of calculations 
is very simple, namely 

1 ‘M 

- [( AL(d)AL(d)A + mA 2 L(d)A 2 L(d)A 2 ),L(d)] = - ^ □ [w 2 , /^]<9A 

Most of the terms in expression (j2.14j) in calculating the commutator vanish. The only term of 
the third order in partial derivatives survives by virtue of [co 2 ,^] ^ 0. 

Finally, we note that one can get rid of the matrix u ; 2 before the higher order derivative if 
instead of the initial equation (12.11) (more exactly, its left-hand side) one considers the most 
general and more symmetric expression 

~{(AL(d)AL(d)A + mA 2 L(d)A 2 L(d)A 2 ),L(d)} (2.16) 

= ^ - \ ({ft,,M + - {ft,, Mu, 2 )d^a- - m 2 I. 

The matrix uj 2 in the first term disappears by virtue of the property (1A.2D . We cannot eliminate 
this term within the framework of the Duffin-Kemmer-Petiau formalism in principle. On the 
other hand its structure up to a numerical factor coincides with the corresponding term on the 
right-hand side of the equation (II.140 . Here the other question arises whether one could remove 
the term of the second order in cP in Eq. (12.161) . The remainder of the paper will be devoted 
to answering this question. 
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3 Cubic root of the third order wave equation 


Let us consider now the construction of the cubic root of the third order wave equation. It 
is clear that the “naive” representation of the cubic root as was defined on the left-hand side 
of expression (11.1411 is unsuitable. Even with the use of the most general representation (the 
left-hand side of (I2.16|l ) the undesirable term of the second order in the derivatives survives. 
Besides, the coefficient of the operator 

- 

m 

differs from the corresponding coefficient on the right-hand side of Eq. (11. 14ft by the factor 
(—1/8) and to correct it is also by no means easy. This imply that we cannot get rid of the 
unwanted term and correct the coefficient mentioned above by making use of the properties 
of the matrices A and /3 M only. Here, it is necessary to involve some additional considerations 
of algebraic character. In this section we attempt to outline a general approach to the stated 
problem. 

Let us introduce the following deformed commutator, 

S^ = A^-zP„A=[A,M„ (3.1) 


where z is an arbitrary complex number and perform an analysis of the following expression 


- Am) 3 = -i(E ( l f ) Zl z) z[ z) )d^d u d x - A 3 m 3 

+ m(sWs?U + A~f 3« + sjfUsW)^ 

+ im 2 (Ej /M 2 + AE^A + A 2 ~S z) )d^. 

First, we consider the contribution linear in d M in (13. 2D . namely 

im 2 (E&A 2 + A E+ A 2 Sjf>) & 1 , 

where by virtue of the definition (13.111 we have 

3<;>A 2 = A/}„A 2 - z - 
^ m 

3l 2 =w = -^-zA%A, 

AZ^A = A 2 p,A-zAf3 fl A 2 . 

A sum of these three expressions gives 

(■ Afi^A 2 + A 2 f3^A H- fiff) — z [Afd^A 2 + A 2 /3^A 4-/^). 

v m' y m ' 

We see that for any value of the parameter 2 this expression vanishes by virtue of (12.411 . 
Let us consider the contribution quadratic in d p : 


(3.2) 


m(S <f>3<f)/l + 31SW=W + EfAE^ST. 


(3.3) 
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Note that we have written out the expression (13.31) with no explicit symmetrization with respect 
to the vector indices /i and v. Analysis of the expression in parentheses in (13.31) is now more 
cumbersome. Our first step is to write out explicitly each term in the expression (13.31) . 


= Ap„Ap v A - z^A 2 /3 u A - zAp^A 2 + z 2 ^A^A 2 , 
AZ^Z^ = A 2 ^Ap u - zA^A 2 ^ - zA 2 (3^ v A + z 2 A^A^A, 

Z&A ~ « = A^A 2 ^ - zA^A^A - - zp„p v + z 2 /3 ll A 2 /3„A. 

h- f m 

A sum of these three expressions after collecting similar terms is 

(1 — z + z 2 )A(3 fl A(3 u A — u — z^AftnPvA 2 + A" j3^f3 v A^ 

+ [( — z + z 2 )PfjA 2 p v A + z 2 Pf 1 Ap i , A 2 ] + [(1 — z)A/3 f j,A 2 /3 v -\- A 2 /3 fJ ,A[3 v ^. 

Further the use of the identity 




1 

2 


{APfjP v A 2 + A 2 p il p l/ A) + {Ap^Afi p v + P^A^ p v A^ 


enables us to rewrite the sum in a more compact form 


£(z > \ 


— ftA 

m 


{AP^Pu A 2 + A 2 Pfj.p u A) + (A 2 Pfj.APu + Pfj.Ap u A 2 ) 


Here, we have introduced the function 


e(z) = l + z + z 2 = (z- q)(z - q 2 ), 


(3.4) 


which is of great importance for the subsequent discussion. From the expression obtained we 
see that the quadratic contribution (13.31) may vanish if as the parameter z one takes a primitive 
root of equation z 3 — 1 = 0, i.e., q or q 2 , Eq. (11.71) . In addition, it should be noted especially 
that the expression (13.31) vanishes without any symmetrization over the vector indices. 

Now we need to analyze the term cubic in d M in (13.21) . The initial expression is 

-isWspsi'Wa*, ( 3 . 5 ) 

„( Z ) 

where by virtue of the def ini tion of the matrix Z fl ; we have 

S (d 5 W S W = APfjAp u Ap x - PfjAp v Ap x A - zAPf i p i /A 2 p x - zPfjA 2 p u Ap x (3.6) 


+ z 2 PfjAPvA'* p x — zAPfjAp,/P x A + z 2 APnP u Ap x A + z 2 P^A 2 p v p x A. 

A somewhat lengthy computation has shown (see Appendix [Cj) that the contribution cubic in 
<9Eq. (1 3. 50 . in the choice z = q and symmetrization over the vector indices turns to zero. Here 
we only note that the following two equalities, 


APfjA — {fin ^V^3 £ m ) , 

A 2 fifjA = - — {Pfj + iV3£ffi) 


(3.7) 
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are rather useful in the analysis given in Appendix [Cl In the equalities (13.711 we have introduced 
the matrices setting by definition 

_ 2 

These matrices have already been considered in the paper by Azimov and Ryndin m, where 
they played a role of a spin pseudovector for a spin-1 particle. Besides, the matrices were 
intensively used by Fushchich et al. [28] for establishment of the complementary (non-Lie) 
symmetry of the Duffin-Kemmer-Petiau equation. 


4 The r ] fJL matrices 

Let us analyze the results of the previous section from a slightly different point of view. For 
this purpose we introduce a new set of matrices r/ /t that would satisfy the following condition, 

Aq^ = w q^A, (4.1) 

and as an immediate consequence 

A 2 r 7 ^ = w 2 q^A 2 , (4.2) 

where w is some complex number. We return to the expression (11.141) . Here, on the left-hand 
side, instead of the original matrices we set q^: 

[A^q^d^ — m /)] 3 = —i( K Aq ll Aq u Aq\)d^d u d x — m 3 A 3 

+ m (Aq^Aq v A + Aq,j,A 2 q u + A 2 q^Aq u ) d^d u (4.3) 

+ im 2 (A 3 77 Ai + Aq^A 2 + A 2 q fl A')d fl . 

We use the rules of the rearrangements (14.11) and (14.21) to bring the matrix coefficients preceding 
the partial derivatives into a simple form: 

o 1 

Aq^Aq u Aq x = w — q^q„q\, 
m 

Aq^Aq^A + Aq^A 2 ^ + A 2 q fl Aq v = we(w) — q^q u , (4.4) 

m 

X 

A 3 q fl + Aq^A 2 + A 2 q fl A = e(w) — 7 / M , 

m 

where the function e(w) is dehned by the expression (13.41) . It is evident that if we set the 
complex number w equal to q (or q 2 ), then (14.31) reduces to 

— i — quriv'qx d' J d v d x — m 2 I. (4.5) 

m 

Further, if the matrices q M satished the identity of the form (11.131) we could reproduce the 
right-hand side of the relation (jl . 14j) (with the replacement [3^ by q M ). 

Let us now turn to the construction of an explicit form of the matrices q To this end, we 
return to the generalized commutator (13.11) in which for definiteness we set z = q. We rearrange 
the matrix A to the left 


[A,Pp] q = A^ - q/3^A = A(/3 m - mqA 2 (3 ll A). 
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Here, we have taken into account that A 1 = mA 2 * . On the other hand we can rearrange the 
same matrix to the right 


[A,/3„\q = - q^A = [mA^A 2 - q^)A. 

Finally, with the use of an explicit form of the matrices A/^A 2 and A 2 /3 li A, Eq. (13.71) . we derive 
the final form of two equivalent representations of the g-commutator, 


[A, f3fj,] q — A 


1 + 2 9 ) + 


iy /3 




+ Q ) At + 


iy /,3 


A. 


The expressions in square brackets are related with each other by a simple relation 


+ Q )/3fj, + 


iy /3 


= Q 


1+2 9 )Pij j + 


iy /3 




It is clear that as the matrix r/ /t in (14.11) it is necessary to take the following expressionEI 


Vn= ( !+ 2 q )^ + 


i^3 




(4.6) 


and the complex parameter w should be set equal to q 2 . Thus, the rules of the rearrangements 
of the matrices A and r/ M can be written in the final form, 

Arj^ = Q 2 V^A a2 Vv = QV^A 2 - ( 4 - 7 ) 


In the choice w = q 2 according to (14.41) the linear and quadratic in d^ contributions in (14.31) 
vanish. However, as we already know from the results of the previous section, the contribution 
in (14.5)1 cubic in the derivatives after symmetrization with respect to the vector indices also 
vanishes. Here we can trace in more detail the reason of this strange fact. By using an explicit 
form of the //-matrices, it is not difficult to see that now instead of the identity (II . 13)1 we have 




1 + 5« 





(4.8) 


The expression in square brackets is formally equal to 


1 + 5® 



1 + q + q 2 = lim e{z) = 0. 

z^rq 


(4.9) 


Let the function e(z) be a small but finite quantity. It is clear that the quantity is defined 
correctly to an arbitrary numeric factor. For example, the left-hand side of (14.91) can be formally 
represented as a product of two multipliers 


1 + 2®l - 



1+ 2 q 



1 + \ q 


+ 



1 


2 The notation r we have introduced for the matrices (14.611 . is not quite appropriate. In the general theory 

of the DKP algebra PH2HI3D] usually by this symbol the specific expression, namely = 2/3 2 — g /i/i , is meant. 

However, by virtue of the fact that we do not use these matrices in the text, this should not mislead. The only 
exception is Section 5, where we will need a particular value of the expression written out just above for g, = 0. 
To avoid confusion, we set off the symbol rj 0 (= 2/3q — 1) in bold. 
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= (1 — g 2 )(l — g)(l + q + q 2 ) = 31im e{z) = 0. 

z^-q 

Here, we have used the definition of cubic roots of unity (11.71) . The expression obtained differs 
from (14.91) by the factor 3. 

Instead of the operator ( irj^d M — ml), we introduce the following operator: 

ykz)^ {z)a " - mI ) 

and correspondingly instead of the expression on the left-hand side of (14.31) we put 

where we have introduced the notation 









In the limit z —> q according to the formulas (14 .4 p . the contribution that is linear in behaves 
as e 2 ^{z) —>■ 0 and the one that is quadratic in behaves as e 1 / 3 ^) —> 0. On the strength of 
Eqs. (14.81) and (14.91) . nonvanishing contribution gives us only the term cubic in d M and thus we 
finally obtain the desired expression, 


lim 

z^-q 


A 




U ii^ - 



where 

lim rj^z) = rj^q) = ^ 

Z^rq 

and 7) [l is defined by the expression (14.61) . 


(4.10) 


5 Properties of the ^-matrices 

Let us derive a number of relations to which the matrices satisfy. Our first step is to consider 
the commutator of two //-matrices. In view of the original definition (14.61) we have 

fa, Vv\ = ( 1 + \ q ) &] + y 2 ^ 

+ C1 ( 1 + 2 ^ ^ + ^ ^' 

We recall that = [a/, /3 ^]. By making use of the formulas of the algebra in Appendix [Aj 

it is not difficult to obtain the following relations: 
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and 

whereupon 

i[Vp,Vv] = 
Here, we have denoted 


[f/1, Pv\ + [P», 6] = [<*>, [P/M Pu]], 


1 + 2 « 


ii/3\ 


)« 2 + + 


iV 3 


ks‘«]. 


(5.1) 

In view of the property (1A.8P the last term on the right-hand side of the above expression is 
equal to zero. In the first term by the coefficient preceding the matrix SjfJ, the expression 
(14.911 is meant. Thus, as the commutation relation for the r/-matrices we take the following 
expression: 


lim 


= sjff. 


(5.2) 


( \ v 1 IH' \ / 5 I 1 ' \ / J ^ LLV 

e{z) * 

This relation will be deeply used in the next section in analysis of the interacting case. 

In addition, the relation (15.211 enables us to clear up a question about the relativistic invari¬ 
ance of equation 

i4 (^^( Z )^‘ mJ )^ X;Z ) =0 ( 5 - 3 ) 
(in the notation of the wave function ifj we have explicitly separated out the dependence on the 
deformation parameter z). In fact, let us consider the double commutation relation with the 
rj- matrices. By using (15.2D we have 


lim “TT 

2 -><? £{Z) 




(5.4) 


On the strength of the definition of the //-matrices we write out the right-hand side 

ls%},v a] = (l + \s\S,Dx] + is f fj,u 

The hrst commutator, by virtue of the trilinear algebra of /3-matrices, equals 

[s$M = -Pv9\»), 

and the second one by using the same algebra and the property (1A.8D does 

a ] = v - £vg\n)- 

Gathering the expressions obtained, we finally ford 

[S^M =i{g^g\v - Vugx^)- 

If by analogy with (15.ip we introduce a new matrix S$ ( z ) setting by definition 


(5.5) 


S®{z) = 


e 2 / 3 (; 


■i[Vp(z),Vv(z)\, 


18 

















then the double commutator can be presented in a more customary form 


lim [S${z), 


z^q 


£ 1 / 3 (x 


■V\(z)\ =i(Vn9\v-Vv9\n)- 


This relationship ensures us the invariance of equation for the wave function z ) under a 
Lorentz transformation when we pass to the limit as z tends to q. 

Further, let us consider the question of a trilinear relation to which the matrices have to 
satisfy. In other words, what is analog of the relation (1 1.311 for these matrices? Here, we skip 
the calculational details and give only the final result 


lim ~TT 

e[z) 


(%( z )’l\( z )*lv( z ) + Vv( z )Vy( z )Vp( z )) = 9l*Vv + S-A’Jc- 


(5.6) 


The proof of the trilinear relation is presented in Appendix [Dl 

One more interesting question is connected with the behavior of the 77 -matrices under the 
operation of Hermitian conjugation (denoted by the symbol f). First of all we note that instead 
of the expression (14.1011 for the cube of the first order differential operator, an equivalent 
expression could be used 


lim 

z^q 


m^A 2 (- 


■ — ml 


— ( —i — Orjnd^ — m 2 I 
m 


\ eV\z) 

where the matrices fj M (= lining satisfy the (/-commutation relations 

A 2 Vv = A% = q 2 fj^A, 

and their explicit form is defined by the following expression 


% = (1 + 




The matrices and are related by the simple relationship 

1 -q 




9 


2 % = -ev 


(5.7) 


By virtue of the initial definition (14.6p we have 


% 


1 = i + 5«* d 


■f) 


g’d- 


(5.8) 


On the strength of the properties of cubic roots of unity the equality q* = q 2 holds. Further, 
for the /3-matrices the following relation, 


Pi = VoPnVo 

is true. Here rj 0 = 2(/3 0 ) 2 — 1 (see footnote [2] in the preceding section). Besides, by making use 
of the Hermitian character of the uj- matrix, we find 

- w(r/o^T/ 0 ). 
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Under these circumstances, multiplying the expression (15.8)1 on both sides by r] 0 and taking 
into account the properties 


Vo A = -PiVo, VoPo = PoVo, Vo^Vo = -w, Vo = !, 


we finally obtain 


»7o^»7o = (i + ^ - (^) = V 

Comparing the expression above with (15.711 . we derive the desired rule of Hermitian conjugation 

VovlVo = ~QVn- 

In closing this section we discuss the question of the existence of such a nonsingular trans¬ 
formation T that would connect the matrices P^ with the matrices r/ /t , i.e., 


Tft.Tr 1 = v 


(5.9) 


Let us seek the matrices T and T 1 in the form of an expansion in powers of ui 

r T = a-\-bu-\-coj , 

T — a-\-boj-\-cuj , 

where a, b, ... are some unknown constants. Substituting these expansions into the left-hand 
side of Eq. (15.911 and using the formulas (1 A. 1 )l — (1 A.3 (1 . we derive the first system of algebraic 
equations for unknown quantities 


aa + ac= ^1 + - q ), 
ca — ac = 0, 


ba = q 


ab = —q 


iV 3 


iV 3 


We define a solution of this system as a function of two arbitrary quantities a and c 


a = 


c = 


(a + c) 

c 1 


1 + 2 «>' 


1 + 5«t 


, 1 fiV3 

* = 4rhr 

a \ 2 

q 


b = 


a (a + c) 

Further, we require the fulfilment of the relation 

IT ” 1 = / 

that gives us the second algebraic system: 


1 + 2 « 


(a + c) 


W 3 


(5.10) 


(5.11) 


aa — 1, 

(ab + ba) + (be + cb) = 0, 
(ac + ca) + (bb + cc) = 0. 


(5.12) 
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The use of the first equation in (15.1211 enables us to express all coefficients through an arbitrary 
constant a, 


1 - 1 A 

a a \ 2 / 

, _{iV 3 N 

b — a | —— I q, 


1 1 

c = x r q 
2 a 


c=-a q 


and ipso facto the required transformation has the following structure: 


T = 


T~ 1 = a 


/ + 
I - 


iV 3 


2 

iyj 3\ 


qu + -qu 


jquj + -quj 


However, a straightforward multiplication of these two expressions leads to 

IT -1 = 


(5.13) 


instead of the desired one (15.111) . This points to the fact that there is a contradiction in two 
remaining equations of the system (15.121) . Substitution of the solution (15. 13j) into the second 
equation of (15.121) results in the identity, and the third equation gives 



q + q 2 


- 1 . 


The equation does not vanish. This tells us that there is no a nonsingular similarity transfor¬ 
mation connecting the matrices [3^ with r/ /t and in this sense, they are nonequivalent However, 
it can be supposed that these matrices are related in a somewhat weak (limiting) sense. 


6 Interacting case 


In the interaction free case we have derived the expression for the third order wave operator as 
a limit of cube of a certain first order operator, namely, 


lim 


A 





( 6 . 1 ) 


It remains to take up the question of a modification of this expression in the presence of an 
external electromagnetic field. We introduce the interaction via the minimal substitution: 


d* = d^ +i e A»(x). 

With an external gauge field in the system the left-hand side of Eq. (16.11) takes the form 

r / \ -i 3 


lim 


A 


£ 1 /3(x 


rjn(z) D M — ml 


From the last two terms in (11.31) and from the corresponding relations in (14.41) it is not difficult 
to see that the contributions linear and quadratic in the derivatives vanish in the limit z —)• q 
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in the interacting case also. This is independent of the eventual noncommutativity of D' s and 
in doing so, as in the interaction free case, we have the following expression!!: 

i 3 lirn - (ArinAr] u A7]xD fJ 'D 1 'D x ^\ — A 3 m 3 = — i — lim - ( rinii^rix D u D^] — m 2 . (6.2) 

z^rq £\ J 771 z^q £\ J 

However, here one can already expect that by virtue of noncommutativity of the covariant 
derivative this limit will have overwhelmingly more complicated structure in comparison with 
the right-hand side of (16.1)) . 

For analysis of the expression (16.21) we make use the following identity for a product of three 
covariant derivatives 

6 D ti D u D x = {D», D\ D x } (6.3) 

+ [3 D^{ieF vX ) + D x (ieF^) + D u (ieF^ x )\ + [2 {ieF< iX )D v + 2 (ieF^)D x ], 

where for brevity by the symbol {D 11 , D u , D x } we mean a product of three .D-operators com¬ 
pletely symmetrized over the vector indices /i, v and A: 

{D^,D U ,D X } = (D^D u D x +D x D u D ll ) + (D u D^D x +D x D^D u ) + (D^D x D u +D u D x D^). (6.4) 

The Abelian strength tensor F lw {x) is defined by 

[D< M ,D V ] = ieF^(x). 

The proof of the identity (16.3|) is given in Appendix lEl 

Our first step is to consider the contribution in (16.21) due to the symmetrized part (16.41) . In 
view of a total symmetry over permutation of the indices, we get 

-i — lim - 7) fl 7) u 7)x {£> M , D v , D x } 

m z—>q £ 

= 6m ^ \ + VxVvVn) + (v^Vx + VxVuVv) + (VvVxVn + VuVxV^l { D ^ D v , D x } 

= ~ ll tn ( 9ftl,r,x + 9vXr} » + 9i^xVu)i D ^ D u , D x }. 

In the last line we have taken into account the trilinear relation (15.61) . By making use of an 
explicit form of the symmetrized expression (16.41) and collecting similar terms we obtain finally 
the desired limit 

-i- lim - 7 hVuV x{D^,D u ,D x } = -t - [D 2 ( Vx D x ) + D^ V xD x )D^ + ( Vx D x )D 2 ]. 

m z^q £ m 

Let us next consider the contribution from the expression in the first square brackets in 
(16.31) . The initial expression for our analysis is 

- i — lim - r )ll r) I/ r)x\3D ll (ieF uX ) + D x (ieF^) + D v (ieF^ x )} (6.5) 

771 z^tq £ 

= e TT- lim -^ D ^F uX {37^[r] p ,i]x} + [VvAlxlVv + (VuV^Vx ~ VxV^Vu)}- 

ZjTDj z ^ Q £ 

3 In the subsequent discussion for simplicity the ^-dependence of the various quantities under the limit sign 

is understood although not written out explicitly. 
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For the first two terms in braces by virtue of the property (15.21) we can set at once 

lim -i[g^g u \ = SjS) {=i[P„,p v ]). (6.6) 

z->q £ 

For the contribution in parentheses we use the identity 




(6.7) 


and then the relation (16.61) . Eventually the limit of expression (16.51) takes the form 
~ DUieF^) (4+ S^n„) + F D^ieF-^SfJn, - S^V). 

The expression in the second square brackets in (16.31) is analyzed in just the same way with 
the use of Eqs. (16.61) and (16.71) . Collecting all of the above calculated in (16.21) and recalling the 
factor 6 on the left-hand side of the identity (16.31) . we derive the desired expression for the cube 
of the linear operator in the presence of an electromagnetic field, 


lim 

z—>q 


A 




-i A {2 \l)'‘(ri > J) i ) + l)‘ , ( n> l)'ll) l ,-t n> l/ i )l/‘\ (6.8) 


+ i e [D“F->- (4 + S { fX) - D“F- X (Sj£ Vx - Sf^)] 

+ e[F- x D'‘( % S^ + S^%) -rt'fSW,, -mV. 

The above expression has been presented in the most symmetric form. However, it can be 
rewritten in a slightly different form. This will enable us, in particular, to compare it with 
a similar expression suggested earlier by Nowakowski [I], At the beginning we consider the 
expression in the first square brackets on the right-hand side of (16.81) . In the second term there 
we rearrange the operator D M to the right, 


D' j (t Ix D x )D II: = ( Vx D x )D 2 + r lx [D'\ D x ] 

= (■ )]\D X )D 2 + ^ieF uX D ll {g^ u r] X - g^ x g v ). 

Similarly, in the first term at the same place we rearrange the operator D 2 also to the right 

D\g x D x ) = (rj x D x )D 2 + ^ieF uX D^(g^g x - g^xVu) + ^ieD^F uX (g^g x - g^xVu ). 

Thus, instead of the expression in the first square brackets now we have 

3 (g x D x )D 2 + ieF vX D^(g^g x - g^ x g u ) + ^ ieD^F uX [g^g x - g^ x g v ). 

Further we rearrange the matrix in the terms S^g^ (the second and third square brackets 
in (16.81) 1 to the left 

sUxVi* = Vn S lx + [SlxiVu.] = v»slx ~ g^xVu)- 
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Here we have used the property (15.51) . Gathering all of the above calculated and collecting 
similar terms, finally, we derive instead of (16.81) the following expression: 

r / \ -I 3 


lim 

z^tq 




= -i — (VnD^D 2 -m 2 I 
m 


5 / le 

6 

4 


6 V 2 m J G 5 -)+ 3 £ ) W* ~ D * F 


v\ 


6V^ 


+ J ) (S^Jn, - g^F^D*. 


(6.9) 


Let us compare the right-hand side of this expression with that of Nowakowski (Eq. (15.71) in [lj). 
In our notation, we have here 

- i — (finD^)D 2 - m 2 I 
m 

-(£)GkI , )W' A (6.10) 

(£) (P^F^D" + (W) (SjSh - g^F^DF 


In rewriting this expression we have used the identity for the /3-matrices 

Wl - foPM = - (s'jSth - sffpj. 


The first difference between (16.91) and (16.101) is that in (16.91) in all terms instead of the matrices 
fin we have the matrices r ] M (except for the spin structure = i[fi v , fi\\ 1 which is the same in 
both cases). Besides, in (16.91) the third term has no analog at all. It is also interesting to note 
that the numerical coefficients in the last four terms in (16.91) were generated in such a manner 
that if the operator commuted with the function F uX (x) the structure of two expressions 
(16.91) and (16.101) (up to the corresponding matrices) might have perfectly coincided. 

Thus, in this section we have shown that our approach correctly reproduces the general 
structure of the third order wave operator in the interacting DKP theory as it was suggested 
by Nowakowski in [I]J on the basis of purely heuristic considerations. The expression obtained 
has a more symmetric form in comparison with Nowakowski’s one, Eq. (16.101) . However, the 
main difference of our expression (16.91) from the corresponding expression (16.101) is that here 
instead of the original /3-matrices we need to introduce more complicated combinations, namely 
//-matrices. The physical meaning of this fact so far is not clear for us. 


7 The general structure of a solution of the first order dif¬ 
ferential equation ( 15.31 ) 


In this section we analyze the general structure of a solution of the equation of the first order 
in the derivatives, Eq. (15.31) . With an external electromagnetic field in the system, the equation 
takes the form 


C(z , D)ip(x] z) = 0. 


(7.1) 
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Here, we have introduced a short-hand notation for the first order differential operator 


£(z,D) = A^ y 1 (7.2) 

In the notation of this operator we have explicitly separated out the dependence on the defor¬ 
mation parameter z. A solution of Eq. (17.11) can be unambiguously presented in the following 
form: 

^{x;z) = [£{z, D)] 2 (p(x;z), (7.3) 

where in turn the function (p(x; z ) is a solution of the third order wave equation 

[£(*,£>)]* <p(x’,z) = 0. (7.4) 


Such a representation of the general solution of the main first order equation (17.11) is the most 
convenient in practice by virtue of the fact that the wave function z ) is a regular function 
of the parameter 0 in the limit z —)■ q, whereas ^>{x\z), generally speaking, is not regular (see 
below). The regularity of the function tp(x] z) is a consequence of the existence of a well-defined 
limit of the cube of the operator £: 

lim I" C(z, D)] 3 = rhsof Eq. (16.91) . 

z^rq L J 

Let us analyze in more detail the structure of the solution i()(x;z) in the form (17.31) . For 
simplicity we restrict our attention to the interaction free case. We introduce the following 
notation: 

S = z — q. 

Then one can present the matrice^l r/^z) in the form of an expansion in terms of <5: 

V,m{z) = [^ + \z\^ + z{^Y 4 S ji li = 'n lx + 8rf ll . (7.5) 

Here, the matrices r/ /t are defined by the formula (14.61) . and the matrices r)'^ have the form 

= ( 7 - 6 ) 


Taking into account the expansion (17.51) and the definition of the function e(z), Eq. (13.41) , the 
first order differential operator C(z, d ) can be rewritten as follows: 


£{z,d) = 


A 


^77 -jr /3 ^ dM + ^ 2/3 ^ %d 11 - ml 


g = q-q- 


(7.7) 


From an explicit form of this operator it is clear that a solution of the equation (17.11) can be 
obtained in the form of a formal series in positive and negative powers of the parameter h 1 / 3 : 


^{x; z) = ... + i *0—2/3 (x) + ^-1/3 (®) + ^o( X ) + h 1/3, 0i/3 (x) + ... . (7.8) 

4 It should be noted that the matrices 77 ^( 2 ), generally speaking, are defined up to an arbitrary matrix 

function 7Jt(z) such that lim 2 _> g TA Z ) = 0- However, we set for simplicity 'T tl {z) — 0 throughout this paper. 
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It is naturally to be expected that the wave function 0(x; z) will be singular in the limit z —>• q. 
It can be seen more precisely from analysis of the representation (17.31) . Really, an expansion of 
the square of the operator £(z, d ) has the form 

[£(z,c 9)] 2 = (7.9) 

- ^73 -^73 d " dv ~ -^3 -JTz m ( A ^ A + A V) ^ + m2A - 2 

- £ 1/3 {Av» A vl + Av'^AVv) - 5 2/3 m (Arf^A + A 2 r ] ' fl ) 

Q Q 

- 61,3 ~^J3 (, AT i A ’&) a ‘‘&' 

Further, by virtue of the fact that the solution (p(x; z) is regular at z = q, it can be presented 
in the form of a formal series in positive powers of 5 1 / 3 : 

ip(x; z) = ip 0 (x) + 5 1/3 ^i/ 3 (ic) + £ 2/3 W 2 /z{x) + <Vi(x) + • • • • (7.10) 


Substituting the expansions (17. 8ft - (17.101) into the relation (17.31) and collecting terms of the 
same power in 5 1 / 3 , we obtain that = ^-a/ 3 {x) = • • • — 0 and 

= - {Ar]^Ari u )d lx d l '(p 0 (x), (7.11) 

Q 

% 1 

-0-1/3 (x) =- Y/ 3 m (Ari fl A + A 2 r] IJ )d^ipo(x) - —^ (Arj^Arj,,)d IJ, d u (p 1 / 3 (x), 

X 1 

0o(^) = m 2 A 2 (f 0 (x)- y^m(A'q^A + A 2 7i fJ )d ti <f 1/3 (x) - {Ar]^A^ u )d^d l/ ^ 2 / 3 {x) } 

and so on. Thus, if <-Po{x) ^ 0 and/or <£>1/3(2;) 7^ 0, then the solution of the first order differential 
equation (17.11) is singular with respect to the parameter S 17,3 in the 5—^0 limit. The maximal 
power of the singularity is equal to 2. 

The differential equations to which the functions <po(x), (pi /3 (x), ... must satisfy are defined 
by the corresponding expansion of the cube of the operator £(z,d). With allowance for the 
expressions (17.7[) and (j7.9|) . we get the following: 

1. the singular contributions'- 

r 1 : -~(A Vll Ar lu Ar ]x )d^d u d\ 

5 ~ 2/3 : m (Arj^A 2 ^ + A 2 rj tl Arj u + Arj^Arj^A) d' J d v , 

5 ~ l/3 : — l T j^rn 2 (A 3 r) IJ , +Arj^A 2 + A 2 r) lx A)d fX . 

The hrst expression vanishes by virtue of the nilpotency property: (77 • <9) 3 = 0. The other two 
vanish on the strength of the properties (14.71) and (11.121) . Further, we have 


2 . 


the regular contributions: 

5 ° : - — [ArjpArivArix + Arj^ArjlArjx + Arj^A^Ar]^] d ,J d v d x - m 2 I = U 0 (d), 

Q 

5 1/3 : -^73 rn [{Arj^A 2 r) u + A^AriJ) + (A^A 2 ^ + A 2 r]^Ar]' v ) 

Q 

+ (Arj^ArjuA + Arj^A^A)] d ,l d v = ti 1/3 (d), 


6 2/3 
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Substituting this expansion of the cube of the operator C(z‘,d) and the expansion (17.101) in 
(im we obtain the desired equations for the functions <£ 1/3 0*0 > ■ ■ ■ , 

5° : U 0 (d)<p Q {x) = 0, (7.12) 

<5 1//3 : Mo{d)ipi/3(x) + Ui/3(d)(po(x) = 0, (7-13) 

5 2 ' 3 : .... 


By using the following 


1 A ^ A2 = + 

(v'^lxVu + VvVWv) + {VuVxVv + VvrfxVn) + (VuVWv + rfvVxVn) = q(9^ aW + 9u \Vn) 

and the properties (14.711 . it is not difficult to show that the operator U 0 (d) is reduced to our 
third order wave operator 

Uo(d) — (—i — Orjud^ — m 2 l\ 

\m ) 

The corresponding first “correction” to the operator can result to the following simple form: 


t/i/ 3 (<9) = -p 1/3 [n,,:n v - qv'^ + q 2 vM d^d". 


By this means in this section we have presented a simple scheme of calculating the wave 
function ^(x; z) satisfying the basic first order matrix equation (15.311 . The scheme is based on 
using the solution <p(x] z ) of well-defined third order wave equation (17.411 . We have shown that 
the required solution x] z ) exhibits a singular character in the limit z —>■ q. This singularity 
has a finite (the second) order in the small expansion parameter 5 1 / 3 . The crucial equation in 
all schemes of calculations is the equation (1 7.12 [1 . ft is the solution ipo(x) of the third order 
wave equation that enables us to define a complete solution <p(x] z) by means of the relations 
of (17.1311 type and then via the relations of (17.1 lH type to define a complete solution ip{x\ z) 
with any degree of accuracy in the parameter h 1 / 3 . The generalization of the results of this 
section for the case of the presence of an external electromagnetic field in the system under 
consideration is straightforward. 


8 The Fock-Schwinger proper-time representation 


In this section we discuss in more detail another difficulty (it has already been mentioned to 
some extent in Introduction) closely related to noncommutativity of the Duffin-Kennner-Petiau 
operator in the presence of an external electromagnetic held 

Ldkp(D) = — ml , (8.1) 

5 The first pair of the relations is a direct consequence of (13.71) and (Id. 71) . The third relation is easiest to obtain 
from Eq. (15.61) by differentiating with respect to z and setting then z = q or by a straightforward calculation with 
the use of the original definitions of the r)- and //-matrices (Eqs. (14.fill . (17.61) 1 and of the relations (ID. 21) (ID. 41) 
from Appendix [DJ 
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and the proper divisor dpxp(D), Eq. (11.51) among themselves. This difficulty is associated with 
the impossibility of defining the path integral representation for the spin -1 particle propagator 
interacting with a background gauge field within the standard DKP theory only. To understand 
why this is so, we turn again to the Dirac theory. For the spin-| case there are a number of 
well-developed techniques of deriving the path integral representation for the Green’s function 
of a spinor particle in background Abelian mm or non-Abelian mm gauge fields. Our 
main interest here is with the very first step in such a construction. It is connected with the 
Fock-Schwinger proper-time representation of the inverse Dirac operator L Dirac( D ) = 
ml ) _1 . This step consists in “squaring” the denominator through multiplying the latter by the 
corresponding Klein-Gordon-Fock divisor 

1 _ _ dpi rac (D) _ ^ 2) 

LpiracyD) dp irac (D)Lp irac (D) 

where dpi rac (D) = (iy^ZA + ml) followed by the Fock-Schwinger proper-time representation 
(see below). It is worthy of special emphasis that the basis for this “obvious” passage is a 
simple, but very important fact: commutativity of the operators dpi rac (D) and Lpi rac (D) among 
themselves. 

Let us return to the Duffin-Kemmer-Petiau theory. Given the explicit expressions for the 
operators L dkp (D ) and d DKP (D), Eqs. (18.11) and (11.51) correspondingly, by analogy with the 
Dirac case, we seemingly could write at once 

1 _ dpxp(D) 

Lpkp{D ) dp KP (D)L DKP (D)' 

and further follow the known procedure. However, by virtue of noncommutativity of these 
two operators among themselves, the expression on the right-hand side is clearly meaningless. 
For the construction of the needed path integral representation for the L~j^ kp {D) operator we 
inevitably come to the necessity of introducing into consideration a divisor that would commute 
with Lpxp(D) and eventually result in the third order wave equation. Below we will briefly 
describe our approach to the problem under consideration. 

We return again to the spin-| case. Following the paper by Fradkin and Gitman |9] instead 
of the initial Dirac operator Lpi rac (D) we introduce the operator transformed by the factor 275 , 

C, = C-Piraci^D) Ull). 

We know this operator to be the square root of the Klein-Gordon-Fock operator (more exactly, 
one of its roots), as it is defined by means of Eq. (11.91) with the replacement d p —> D t ,. It 
is precisely for the inverse operator LA 1 that it is natural to determine the Fock-Schwinger 
proper-time representation. We believe the operator C to be an odd (Fermi) one by definition. 
Instead of (18.2)) we have now 

^ OO 

j = £- = iJdrJ^ e ~ir(H - ie) + r X C £ ^ +0? (8 . 3) 

0 

where 

H = t 2 = -(D 2 + m 2 )l + | ± 
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r is an even variable and x is an odd (Grassmann) variable, anticommuting by definition with 
£. The pair (r, x) is treated as a proper supertime. By virtue of the fact that the Hamilton 
operator H is represented by a product of two Fermi operators, it is an effective Bose operator 
as it must be. 

For the case of DKP theory as the operator £ we take the cubic root of the third order 
wave operator, namely the expression (17.21) . Let us assume that this operator is a para-Fermi 
operator (parastatistics of order two). In this case it is not difficult to write an analog of the 
Fock-Schwinger proper-time representation for the inverse operator £ _1 similar to (18.31) . 


1 

I 


£}_ 

£ 3 



-ir(H(z) - in) + i (r[x,£] + jT 2 [x,C] 2 ) 


(8.4) 


where now 

H{z) = £ 3 (z, D) 

and x is a para-Grassmann variable of order p — 2 (i.e., x 3 = 0) with the rules of an integration 

Jd 2 x = 0 = Jd 2 x [x, £], Jd 2 X [x, t] 2 = 4 £ 2 . 

We consider that the para-Grassmann variable x anc i the operator £ conform to the following 
rules of commutation: 

[[x,£],£] = 0, [[x,£],x] = 0. 

As a proper para-supertime here it is necessary to take a triple (r, x, X 2 )- The Hamilton operator 
H{z) represents a product of three para-Fermi operators; therefore, in this case, too, the H(z) 
is an effective Bose operator. Of course, only its limiting value has a physical meaning 


H = lim H{z) 

z—>q 

where the most right-hand side limit is defined by the expression (16.91) . 

However, it should be specially noted that the expression (18.41) for an arbitrary value of the 
deformation parameter 0 is meant here as a purely formal one, since the fact of the presence of 
supersymmetry corresponding to parastatistics of order two has not been demonstrated by us 
explicitly. In addition, the indicated parasupersymmetry most likely does not take place for z 
distinct from q (or from q 2 ). The final conclusion about the existence of this symmetry can be 
made only after constructing the appropriate path integral representation and passage to the 
limit z — y q. In this connection one should point out a similar situation taking place in the 
case of the so-called deformed Heisenberg algebra with reflection R introduced by Wigner (34]. 
The algebra includes a real-valued parameter //, and reveals very peculiar properties for special, 
discrete values of this parameter (y = —(2 p + 1), p = 1, 2, ...). This very nontrivial fact was 
first found by Plyushchay (35], who pointed out the relationship of the (2 p + 1) - dimensional 
representations of the R - deformed Heisenberg algebra to parafermions of order 2 p. These spe¬ 
cial values of the parameter v also reveal themselves in the context of supersymmetry as well 
as field theory [36]. 


= lim £ 3 (z, D) = lim 

z^tq z—^q 
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Besides the presence of (local) parasupersymmetry of order two in the system under con¬ 
sideration will ensure the consistency of the minimal coupling prescription for the case of the 
DKP theory considered in our paper, as it is, for example, for the spin-1/2 charged field. In the 
latter case a local supersymmetric structure associated with the Dirac equation guarantees a 
consistency of the minimal coupling prescription when a free theory is generalized for the case 
of interactions with an external electromagnetic held. 

The expression (18.41) can be taken as the starting one for the construction of the desired 
path integral representation with the use of an appropriate system of coherent states in a close 
analogy with the approach proposed by Borisov and Kulish [8] for the spin- 4 case. Here, it 
would be possible to make good use of the known for a long time |37hl39] connection between 
the trilinear algebra of /3-matrices and the para-Fermi algebra of order two. In fact all appara¬ 
tus needed for the construction of the path integral representation (coherent states, formulas 
of orthonormality and completeness, and so on) can be found in the papers by Kamefuchi and 
coworkes [33]. However, in our case instead of the matrices /3 M we have the matrices r/ /t (^). The 
trilinear relations for these matrices formally coincide only in the limit z —>• q. As a hint of what 
we shall do in this more complicated situation, the unpublished paper by Dunne [40] can serve. 
In the latter it was shown how one can define the creation and annihilation operators explicitly 
depending on the deformation parameter z and the corresponding relations of commutation 
with the subsequent passage to the limit z —>• q resulting in the finite expressions. 

In closing let us mention another fact closely related to the subject matter of this section. 
In the literature there are very few papers dealing with the problem of construction of an action 
for a relativistic classical spinning particle using the para-Grassmann variables with subsequent 
quantization of the classical model [411 - 143] . Here we note only that in the action suggested in 
these papers there are the linear and quadratic in para-Grassmann variable x (in our notations) 
terms, which are similar to those in the exponential function in the expression (18.41) . In our case 
these terms automatically appear in defining the Fock-Schwinger proper-time representation, 
and in the works [4T] HU] they insure the invariance of the action under the local world line 
para-SUSY transformation. However, the kinetic part of the action in [411143] was chosen in a 
complete analogy with the kinetic one for the classical models of a Dirac particle, whereas we 
expect based on the general formula (18.4[) that the situation here may be more complicated. 


9 Conclusion 

In this paper, we have set up the formalism needed to construct a cubic root of the third order 
wave operator within the framework of Duffin-Kemmer-Petiau theory. One of the key points 
in our approach is the introduction into consideration of the so-called deformed relation of 
commutation, Eq. (13.11) . On the basis of the latter a new set of the spin matrices was defined 
instead of the standard DKP matrices (3^. It was shown that the third order wave operator 
is obtained as a formal limit of the cube of a certain first order differential operator. This 
operator is singular with respect to the deformation parameter z when the latter approaches 
the primitive cubic root of unity q. Finally, we suggested a way to apply the derived expression 
for the cubic root to the problem of the construction of the path integral representation for the 
Green’s function of a spin-1 particle in an external electromagnetic held. 
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A few words may be said here about the para-Grassmann variables, which will be used in 
the construction of the desired path integral representation. Although the para-Grassmann 
algebra of order p = 2 is still quite visible for concrete calculations, however, probably in the 
situation under consideration the use of its bilinear version |44ff49] (sometimes it is named the 
generalized Grassmann algebra) is more suitable. It is connected with the fact that on the one 
hand the primitive cubic root q explicitly enters into the definition of the p- matrices, into the 
commutation relations and so on, and on the other hand the use of a primitive nth root of 
unity (in particular for n — 3) is directly laid in the basis of the new para-Grassmann calculus. 

It only remains for us to say a few words about the massless limit of the third order 
wave operator. Throughout this paper we have considered that the parameter m, which by 
convention is responsible for the mass of a particle, was not equal to zero. As it is known [26j 
in the massless variant of DKP theory the scalar “mass matrix” ml must be replaced by a 
singular matrix Mcu 2 , where M is an arbitrary constant with the dimension of mass. The 
matrix Mu 2 does not commute any more with everything, as this occurs for ml (see the text 
after Eq. (11.131) ). The last circumstance qualitatively changes the whole picture of calculations 
in comparison with the massive case. Thus, for example, here there is no analog of the formula 
(11.111) . Therefore, now we can only speculate that in the massless case instead of the expression 
(11.141) there should be something of the type 

[B(ip^ - Mu 2 )] 3 = -i T n^d", 

where B is a certain matrix depending on to and additional parameters, with a possible replace¬ 
ment of the matrices by a more complicated combination. Preliminary consideration has 
shown that, probably, one of the crucial factors here is the use of the cubic roots of minus unity 
rather than of unity. This problem requires separate careful consideration. 
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Appendix A The matrix algebra 


In this Appendix we give some necessary formulas of the 

matrix algebra for the spin-1 

case, which are used throughout in the text. Details of the proof of these formulas and also 
their generalizations to higher dimensions can be found in the papers by Harish-Chandra [26] 

and Fujiwara [29]. We use the metric g /w = diag(l, —1, —1, —1). 
the hJ matrix: 

Let us recall the definition of 



Then in view of the above definition and the trilinear relation for /3-matrices, Eq. (11. 3j) , we have 

to 3 = u, 

(A.l) 

^ 2 A + A^ 2 = A> 

(A.2) 

UJ^UJ = 0, 

(A.3) 

fitifii/Ld ~\r , 

(A.4) 

w 2 AA = AA<A 

(A.5) 

A W A + A^A = o. 

(A.6) 

The next formulas 


{A> A} ^ T A’ A} 

(A.7) 

[A>AA~AA>A] = ° 

(A.8) 


are an obvious consequence of (IA.4j) . If one defines the matrix B = then the following 

relations are also valid: 


uj 2 = 3 — B : Bui = uB = 2 uj. 

Besides the useful contractions are 

^u% = Z{\-u% AAA = A- 


(A.9) 


Appendix B Construction of the matrix A 


We write down once again an explicit form of the matrices (A, A 2 , A 3 ) obtained in section 2 

A = a 


T ^.V 3 3 2 

1 + i — u uJ 

2 2 


A 2 = a 2 


A 3 = a 3 1 = — I. 


T .V3 3 2 

1 — i — u UJ 

2 2 

1 


(B.l) 


m 


An immediate consequence of (IB.ID is the following relation: 


- A + \ A 2 + \ A 3 = 3{I - u 2 ). 

a or a 


(B.2) 
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Let us construct such a matrix A that simultaneously satisfies two requirements: 


1 1 j9 1 - 

— A H— t: A* H—^ — 0 

a cr a- 1 


and 

* 4 3 = -/. 

m 


We search for this matrix in the following form 


A = A + x(I — uj 2 ), 


where x is an unknown parameter. It follows that 

A 2 = A 2 + x(2a + x) (I — oj 2 ) 


(B.3) 

(B.4) 


(B.5) 


and 


„4 3 = — I + x{x 2 + 3 ax + 3 a 2 )(J — to 2 ). 


m 


(B.6) 


Substituting the obtained expressions into the left-hand side of the expression (IB. 3D and con¬ 
sidering (IB. 2D . we obtain 


— A A —— A~ ~l~ — n A 2 — —r (ic 3 -|- 4 ax 6a x -I- 3cr 3 ) (/ — uj ). 

a cr a 6 a 6 

The requirement of vanishing the expression on the right-hand side results in the following 
equation for the parameter x 


x 3 + Aax 2 + 6cc 2 a; + 3a 3 = (x + a)(x 2 + 3ax + 3a 2 ) = 0. 

It is apparent that the given algebraic equation has three roots but only two of them are 
compatible to the additional requirement (1B.4D . Actually, by virtue of (1B.6D we have another 
equation for the unknown parameter x 


x{x 2 + 3 ax + 3a 2 ) = 0. 


The two roots needed are 



Thus, we find the required expressions for a new set of matrixes (A, ^l 2 ,^4 3 ) : 

A = A + x(I — uj 2 ), 

A 2 = A 2 + ax* (I — uj 2 ), 



m 


In deriving the expression for *4. 2 here, we have considered in (IB.5D the identity 

x[2a + x) — ax*. 
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Appendix C The proof of vanishing ( 13.51 ) 


By using the relations 02.31) and 
in 03.61) in an identical form: 


one can rewrite the first, second, fifth and seventh terms 


Ap^Ap^Ap x = -Apftp u A 2 p x - Afi^PvPxA, 

P^Ap u Ap x A = -p^pvPxA - P ll (5 v A 2 l3xA, 

AP^APxA = -A^PvPxA 2 - Ap„p v A 2 Px, 

PnAp u A 2 Px = —P^A 2 p v Ap \ - PfiPvPx- 

m 

Substituting these relations into 03.61) and collecting similar terms, we obtain then 

— — z 2 PuPvftx + (ft/ifivA 2 fixA — AP fl A 2 p i/ Px) — {z + z 2 )(3 /1 A 2 (3 u Af3x 

+ {zA 2 ^MxA - z 2 Ap^ v p x A 2 ) + £(z)(P,A 2 (3 u (3xA - Ap,p u A 2 Px). 

The last term here can be turned into zero if one sets z = q. Taking into account another 
identity 

— A(3 fl A 2 f3 l/ (3x — A 2 P^AfivPx H- P^PuPx 

m 

and the equality q + q 2 = —1, we derive further 

1 

m 

+ {A 2 (3ft A (3 u (3 x + PftA 2 p v APx + P/iPvA 2 PxA '). 

The final step is a contraction of this expression with d^d u d x . Making use again of the relation 
q = — 1 — q 2 , the identity 01.131) and the property 02.41) . we obtain 


(1 - q 2 )/3ft/3 u /3x + {qA 2 Pf i p v PxA - q 2 APftP u PxA 2 ) 


1 


- up^ - □ (A%A)d» 


m 


(C.l) 


+ (A 2 f3fjAf3 u f3x + PftA 2 p v APx + P.P.A 2 p x A)d^d u d x . 

Now we are concerned with an analysis of terms containing the matrices A and A 2 . With 
allowance made for 03.71) . the expression (1C.lj) turns to 


1 UPftd* - aL p i J3J3xd> l d"d x + J- UPftd^ 


m 


(C.2) 


2 ^ iy/3 {i^PvPx + PuivPx + PpPui a ) d^d u d x + ^ 


By using once more the identity 01.131) we see that the first three terms here mutually cancel. 
Furthermore, it is not difficult to verify that the following equality holds 


ZuPvPx + P^uPx + P^PviX = uPftPvPx ~ PftPvPxU. 
In view of PC.31) the last but one term in PC.2D takes the form 

- ^ (uPftPvPx ~ P lx P u Pxu)d IJ 'd I 'd x = - iV 3 

This term is canceled precisely by the last term in PC.2D . 


(C.3) 
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Appendix D Trilinear relation for the ^-matrices 


By a direct multiplication of the matrices r/ M , Eq. (14.61) . we derive the starting expression 

.3 


VuVvV A + VxVvVn = 1 + o 9 } (A PvPx + PxPu A) 


(D.l) 


+ |i+2<! 


1 Wi\/3 


+( 1+ 2 g 


2 

ix/3^ 2 


+ ZtiPvP a) + (/i <h- A)] 
<T [(^m&'/^A + CuPvCx + P^u^x) + (/i o A)] 


i«s/3 


q 3 (€n€v€x + £aAA/J- 


For the first term on the right-hand side we use the basic relation for the /3-matrices, Eq. (11.31) . 
For the second term it is necessary to use the relation (1C.3D , which in view of Eq. (11.31) leads to 


(AA/£A + P/xCuP. A + CuPvPx) + (/i 44 A) — (7 M y^A + dXu^fi- 
In analysis of the third term in (ID.ID we first note that 

Pv}^ 9fiv^ P^Pv 


(D.2) 


and 

then 


£,A£a = 9hvU 2 Px + 9xuP^ 2 ~ {^PuP^Px + P^PxPvU 2 ) - 


(CuCuPx + x + P^CvCx) + (/i 44 A) — —2 {P^PuPx + /3a A A) 


(D.3) 


4-uh 


(A, ^}/?a + {/3a, A} A — Pv{Pn, /3 a } 


+ 


/3a(A, A} + A {/3 a , A,} — {/3 m , /3 a }/3„ 


oj 


2 (g^vPx 4” gxvPp) 4~ tu (g^Px 4“ gxvPp) 4~ (g^Px 4~ gxv Pi_p^~ 
^9nvPx 4~ 9xi/@n)- 

In the last step we have used the property (IA.2D . 

Finally, for the last term in (ID.ID we have 

4 " ^xPuiP^^Px 4 " Px^PP) ^{g^iuPx 4 ~ gxuPp) 

4~ (gin/ftx 4~ gxvPp)^ = (s7«Aa 4~ gxu^p)- 


(D.4) 


Here, for the first term on the right-hand side we have used the property (1A.6D . Gathering the 
expressions calculated above and collecting similar terms, we obtain instead of (ID. ID 


1 + l q 


^V3\ 2 

2 I 


V^VuVx + VxVvVv 

l + ^ qj (W 3 a + gxvP?) + ?(0 fits £a + gxuC.fi) 

= ^{q){g^vx + gxuV„), 


whence it follows the trilinear relation (15.6D . 
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Appendix E Proof of the identity ( 16.31) 


Let us present a product of three covariant derivations D M in an identical form 

D^D V D X = {D^, D v , D x } (E.l) 

- (D X D U D» + D V D^D X + D X D^D V + D^D X D V + D u D x D tl ). 

By the symbol { DD u , D x } one means a completely symmetrized expression defined by the 
formula (16. dl) . We transform each term in parentheses in such a manner so that the expression 
obtained contains the term D^D U D X . Here, we have 

D X D U D^ = D^D U D X -D^[D\D X ] - [D^, D X ]D V — D X [D^ D u ], 

D V D^D X = D^D V D X - [D'\ D U ]D X , 

D x D tx D u = D fX D u D x - [D^D U ]D X - D v [D», D x ], 

D^D x D u = D^D V D X - D^[D U , D x ], 

D V D X D ^ = D^D U D X — [D^D U ]D X - D U [D^,D X }. 

Substituting these expressions into (IE.IK . taking into account the equality 

[D»,D V ] =ieF^(x), 

and collecting similar terms, leads to (16.31) . 
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